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Abstract 



An analysis of traveling wave solutions of partial differential equation (PDE) systems with 
cross-diffusion is presented. The systems under study fall in a general class of the classical 
Keller-Segel models to describe chemotaxis. The analysis is conducted using the theory of the 
phase plane analysis of the corresponding wave systems without a priory restrictions on the 
boundary conditions of the initial PDE. Special attention is paid to families of traveling wave 
solutions. Conditions for existence of front-impulse, impulse- front, and front-front travelling 
wave solutions are formulated. In particular, the simplest mathematical model is presented 
that has an impulse-impulse solution; we also show that a non-isolated singular point in 
the ordinary differential equation (ODE) wave system implies existence of free-boundary 
fronts. The results can be used for construction and analysis of different mathematical 
models describing systems with chemotaxis. 
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1 Introduction 

In this paper we study one-dimensional traveling wave solutions for the models in the form 



Here /i > is a constant; f(U,V) and g(U,V) are functions whose properties will be specified 
later; U = U(x,t), V = V(x,t); in the following we put (jl = 1 without loss of generality. 

The model JT]) is known as a particular case of the classical Keller-Segel models to describe 
chemotaxis, the movement of a population U to a chemical signal V (see, e.g., [HO [3]). In system 
(PQ) fi denotes the constant diffusion coefficient; f(U,V)/U is the chemotactic sensitivity, which 
can be either positive or negative; g(U, V) describes production and degradation of the chemical 
signal; it is customary to include also in the second equation of (P) the diffusion term of the 
form DV XX which would describe diffusion of the chemical signal, but, we adopt hereafter, that, 



U t = (»U x -f(U,V)V x ) 
V t = g(U,V). 
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to a first approximation, D can be taken zero. For biological interpretation of the solutions of 
(PQ) we refer to the cited literature, references therein, and to Section 3.2; macroscopic derivation 
of equation ([T]) can be found in, e.g., [I]. 

The chemotactic models are the partial differential equations (PDEs) with cross-diffusion 
terms; these systems possess special mathematical peculiarities [5]. Such systems were used, 
e.g., to model the movement of traveling bands of E. coli [2j[3j[6], amoeba clustering [7J, insect 
invasion in a forest [H] , species migration ^9] , tumor encapsulation and tissue invasion [1(J|, [TT] . 
Many different spatially non-homogeneous patterns can be observed in chemotactic models, for 
a survey see, e.g., [12] and references therein. One such pattern is that of traveling waves which 
spread through the population. 

A traveling wave is a bounded solution of system (JTJ) having the form 

U(x, t) = U(x + ct) = u(z), V(x, t) = V(x + ct) = v(z), 

where z = x + ct and c is the speed of wave propagation along x-axis; u(z) and v(z) are the 
wave profiles (u-profile and f-profile respectively) of solution (U(x,t), V(x,t)). 
Substituting these traveling wave forms into ([T]) we obtain 

cu' = (v! - f(u,v)v')', 
cv' = g{u,v), 

where primes denote differentiation with respect to z. 

On integrating the first equation in the last system we have the wave system of ([I]): 

u' = cu + f(u,s)g(u, s)/c + a, 
v = g{u,s)/c. 

Here a is the constant of integration that depends on the boundary conditions for U(x,t) 
and V(x,t). In various applications it is usually possible to determine this constant prior to 
analysis of the wave system. For instance, considering system JI| as a model of chemotactic 
movement [2], where the variable U(x,t) plays the role of the population density and V(x,t) 
is an attractant, one usually supposes that u(x, t) dx should be finite, which implies that 
a = (e.g., PQ). On the contrary, in our analysis we do not specify the boundary conditions for 
dl]) and consider a as a new parameter. 

Each traveling wave solution of (JTJ) has its counterpart as a bounded orbit of (J2J for some a; 
in our study we elucidate the following question: for which a there exist traveling wave solutions 
of ([T|) and describe all such solutions. We also note that the case of a = does not exclude 
a model with infinite mass of U(x, t) if the traveling wave solution is a front; moreover, the 
solutions corresponding to finite mass can be only impulses (see below for the terminology). 

It is worth noting that due to specific form of system (pQ) with cross-diffusion terms the 
wave system has the same dimension as the initial system (pQ), which significantly simplifies the 
analysis. This is one of peculiarities which distinguishes cross-diffusion PDEs from those with 
only diffusion terms (see also [51 E3] ) • 

We shall study possible wave profiles of ([I]) and their bifurcations with changes of the pa- 
rameters c and a by the methods of phase plane analysis and bifurcation theory |14t I15j . In 
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this way, the problem of describing all traveling wave solutions of system ([T]) is reduced to the 
analysis of phase curves and bifurcations of solutions of the wave system ([2]) without a priory 
restrictions on boundary conditions for ([T]). 

There exists a known correspondence between the bounded traveling wave solutions of the 
spatial model (pQ) and the orbits u(z), v(z) of the wave system ([2]) (e.g., [U [HI dTJ [18] ) that we 
only list for the cases most important for our exposition. 

Proposition 1. 

i. A wave front in U (or V) component corresponds to a heteroclinic orbit that connects 
singular points of ((2]) with different u (or v) coordinates (FigU]a); 

ii. A wave impulse in U (or V) component corresponds to a heteroclinic orbit that connects 
singular points with identical u (or v) coordinates (Fig\T$)) or to a homoclinic curve of a singular 
point (u, v) of <m (FigW). 



(a) 



(b) 
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Figure 1: Correspondence between bounded traveling wave solutions of system ([I]) (on the right) 
and the phase curves of the wave system ([2]) (on the left); the black dots are singular points of 
(|2|). (a) A front-front solution; (b) a front-impulse solution; (c) an impulse- impulse solution 



Hereinafter we shall adopt the following terminology: we define the type of a traveling wave 
solution of (H]) with a two word definition; e.g., a front-impulse solution means that ^-profile is 
a front, and ^-profile is an impulse (the order of the terms is important). 

For system (pQ) several results on the existence of one-dimensional traveling waves are known; 
see, e.g., [H El El El EH ESI [23 Ell 122] • In most of these references the analysis is conducted 
using a particular model which is given in an explicit form. Quite a different approach was used 
in pQ where the authors consider more general model than ([!]) and do not restrict themselves to 



3 



analyzing a model with specific functions f(U, V), g(U, V); instead their aim was to understand 
how these functions have to be related to each other in order to result in traveling wave patterns 
for U and V. We consider a general class of models as well, and our task is to infer possible 
kinds of wave solutions under given restrictions on f(U, V) and g(U, V). 

Our main goal is as follows: we impose some constrains on the functions f(U,V), g(U,V) 
and study possible traveling wave solutions with increasing complexity of /, g. Special attention 
is paid to the families of traveling wave solutions such that the corresponding wave system 
possesses an infinite number of bounded orbits. We present the simplest possible models in the 
form ((T|) that display traveling wave solutions of a specific kind. 

The main class of the models we deal with is defined in the following way. 

Definition 1. We shall call model (H|) the separable model if 

(CI) f(u,v) = f 1 (u)f 2 (v), g(u,v) = g 1 (u)g 2 (v), 

where fi(u), g\{u) are smooth functions for a ^ u < oo; g2(v) is smooth; f2(v) is a rational 
function: 

for b ^ v < oo; here a,b > — oo are real constant. 

The separable model will be called the reduced separable model if (CI) holds and 

(C2) f2(v)g 2 {v) = const. 

We organize the paper as follows. In Section [2] we present full classification of traveling wave 
solutions of the reduced separable models; we also specify necessary and sufficient conditions for 
these models to possess specific kinds of traveling waves. Section [3] is devoted to the analysis of 
the separable model; we show which types of traveling waves can be expected in addition to the 
types described in Section [2] we also analyze a generalized Keller-Segel model, which does not 
belong to the class of the separable models but display a number of similar properties together 
with essentially new ones. Section 0] contains discussion and conclusions; finally, the details of 
numerical computations are presented in Appendix. 



2 Wave solutions of the reduced separable model 

In this section we present the full classification of possible traveling wave solutions of system 
(PQ) that satisfies (CI), (C2). The reason we start with the reduced separable models is twofold. 
First, there are models in the literature that have this particular form (see, e.g., \20 \ 1221 [23] ); 
second, the special form of the wave system allows the exhaustive investigation of traveling wave 
solutions of ([I]). 

The wave system of the reduced separable model reads 

u' = cu + f 1 (u)g 1 (u)/c + a = h(u), 
v' = gx{u)g 2 {v)/c, 

where the first equation is independent of v. 
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2.1 The cell structure of the phase plane of system (I3J) 

We start with the case of general position. We assume that the following conditions of non- 
degeneracy are fulfilled (later we will relax some of these assumptions): 

(Al) h(u) has no multiple roots; 

(A2) <72 (v) has no multiple roots; 

(A3) gi(u) and h(u) have no common roots. 

Traveling wave solutions of ([T|) correspond to bounded orbits of ([3]) different from singular 
points. Due to the structure of system ([3]) it is impossible to have a homoclinic orbit or a limit 
cycle in the phase plane of ([3]), which yields that it is necessary to have at least two singular 
points of ([3]) and a heteroclinic orbit connecting them (see Fig. [Th 5 b) to prove existence of 
traveling wave solutions of ([I]) satisfying (CI), (C2). 

In general, smooth functions h(u), g(u,v) can be written in the form 

h(u) = h(u)(u — ui) . . . (u — u m ), h(u) ^ for any u, 
g(u,v) = gi(u)g 2 (v)(v - vi) . . . (v - v n ), g 2 (v) / for any v. 

The following proposition holds for neighboring roots of h(u) and g2(v). 
Proposition 2. 

i. Let the wave system ([3]) satisfying (^4l)-(^43) have singular points (u,vi) and (u, V2), where 
i>i and V2 are neighboring roots of g2(v) then one of these points is a saddle and the other one 
is a node. 

ii. Let the wave system (|3|) satisfying (A1)-(A3) have singular points (ui,v) and (u2,v), 
where u\ and U2 are neighboring roots of h(u), then these points can both be saddles, nodes or 
one is a node and another is a saddle. 

Proof. Let (n, v) be a singular point of (|3|). The eigenvalues of this point \\(u,v) = h'(u) and 
\2(u,v) = gi(u)g' 2 (v) are real numbers (henceforth we use prime to denote differentiation when 
it is clear with respect to which variable it is carried out). This implies that singular point (u, v) 
of system (|3|) cannot be a focus or center. 

i. The claim is a simple conjecture of condition (A2). 

ii. Let us consider two equilibrium points (ui,v) and (ii2,v). The eigenvalues Xi(ui,v) and 
Xi(u,2,v) have opposite signs due to (^41). 

Consider another pair of eigenvalues X 2 (ui,v) and \ 2 (u2,v) and assume that (A3) holds. If 
the number of roots of gi(u) located between u\ and u 2 is even (or zero) then the signs of these 
eigenvalues are the same. This implies that one of the equilibrium points is a saddle whereas 
the other one is a node. If the number of roots of g\(u) located between u\ and 112 is odd than 
the signs of these eigenvalues are opposite which implies that both equilibriums are saddles or 
nodes (one node is attracting and another is repelling). □ 

Note, that in case ii of Proposition [2] in order to guarantee that both singular points are 
nodes one should have h! ' (ui)gi(ui)g' 2 (v) > and hi (u 2 )gi(u2)g'2(v) > 0. Due to continuity 
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Figure 2: Two types of orbit cells of system © 

arguments there exists a family of orbits of (0) which tend to one of the nodes when z —* oo 
and to the other node when z — ► — oo. 

Taking into account that straight lines u = Ui, i = l..m and v = Vj, j = l..n consist of orbits 
of system ([3]) we obtain that the phase plane of ([3]) is divided into (m — l)(n — 1) bounded 
rectangular domains whose boundaries are u = &i, i = l..m and v = Vj, j = l..n. We shall call 
these domains the orbit cells. 

Due to Proposition [2] it immediately follows that an orbit cell can be one of the following 
two types (up to 7r-degree rotation) that are presented in Fig. [2j The behavior of the orbits 
inside a cell is completely described by the types of the singular points at the corners of the cell. 
Moreover, any orbit inside a cell corresponds to a bounded traveling wave solution of system 
©• 

Summarizing the previous analysis we obtain the following theorem. 

Theorem 1. The system (JTJ) satisfying (CI), (C2) and (Al)-(A3) possesses traveling wave so- 
lutions 

i. of a front-front type (Fig. [7]a) if and only if the wave system ([3]) has four singular points 
(ui,Vi), i = 1,2, which are the vertexes of a bounded orbit cell and every two neighboring vertexes 
are a node and a saddle (Fig. Wp,); 

ii. of a front-impulse type (Fig. \Tjj) if and only if the wave system ([3j) has two neighboring 
nodes (ui,v) and (u<2,v), (see Fig. 0)). 

Remarks to Theorem [TJ. 

1. In both cases the orbits of system ([3]) that correspond to traveling wave solutions of (H|) 
are dense in the corresponding orbit cell. 

2. System (TjQ) has a traveling wave solution which is a front in w-component and space- 
homogeneous in u-component if and only if the wave system ([3]) has neighboring saddle and 
node with identical n-coordinate, see singular points (tti, v\), (u%, v-z) and (u2, t>i), (^2, ^2) hi 
Fig. [2b. 

It is possible to write down asymptotics for u and v profiles (these asymptotics can be used, 
e.g., as initial conditions for numerical solutions of (H|)). We present these asymptotics only in 
the simplest case. 
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Let us assume that the wave system has the form 

v! = A(u — u\)(u — u 2 ) 



(4) 

v' = —B(v — vx)(v — v 2 )(u — u), 
where A, B > are constant. An explicit solution of is 

u(z) = u 2 + - ^ Tirn* * n -1 ; 

1 + Ci exp {A(m -u 2 )z}' 

B(v 2 -vi) (5) 

. (t>i - 62 (1 + Ci exp {A(ui - u 2 )z}) a 

< Z ) = V 2 + ' 

(1 + Ciex.p{A(ui - u 2 )z}) A + C 2 exp {B(u - ui)(vi - v 2 )z] 

where Ci, C 2 are arbitrary constants. We emphasize here that even with fixed c and a there is 
a two-parameter family of wave profiles. 

Let u\ < u < u 2 and v\ < v 2 . We consider non-trivial profiles {C\C 2 7^ 0). It is straightfor- 
ward to show that u(— 00) = u 2 , u(oo) = and, hence, u-profile is a front; v(— 00) = Ui, -u(oo) = 
«i, and v-profile is an impulse. If u < u\ or u > u 2 then n-profile remains the same and ^-profile 
becomes a front. 

Formulas © can be used as a first approximation for wave profiles of system (pQ) even 
in the case A = A(u),B = B(u,v) where A(u) and B(u,v) do not change the sign when 
u £ (ui,u 2 ), v £ (v±, v 2 ). 

It is worth noting that if u\ = u then 

( \ ~ , vi-v 2 

1 + C/2 

i.e., the boundary of the profile depends on an arbitrary constant; we deal with such solutions 
in the next section. 



2.2 Bifurcations of the travelling wave solutions 

Considering c and a as bifurcation parameters we can relax some of non-degeneracy conditions 
(A1)-(A3). 

First we note that the right-hand side of the second equation of system ([3]) does not depend 
in a non-trivial way on c and a and we will not consider the case when (A2) is violated. In 
general, by varying the bifurcation parameters we can only achieve that either (^41) or (A3) 
do not hold. We shall show that in the latter case new traveling wave solutions can appear in 
system CO) satisfying (C1),(C2). 

First let us assume that (^41) does not hold, i.e., function h(u) has a root u of multiplicity 
m > 1 for some a*,c*, and the wave system ([3]) has a complicated singular point (u,v). The 
system can be written in the form 

u ' = (u-u) m qi{u), 
v ' = ( v - v)g 1 (u)q 2 {v), 

where q\(u)q 2 (v)g\(u) 7^ 0. Then the singular point (n, v) of system ([6]) is either a saddle, a 
node, or a saddle-node [14] . For the first two types of critical points, the structure of the phase 
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Figure 3: The phase plane of system fl5J); m is a root of both h(u) and gi(u). The wave solutions 
corresponding to bounded orbits of the wave system form a free-boundary family 

plane of the wave system was completely described above. In case of a saddle-node the line 
u = u divides the plane such that in one half-plane the singular point is topologically equivalent 
to a node, and in the other one it is topologically equivalent to a saddle; due to the fact that 
u = u consists of solutions of ([3]) this type of singular points does not yield qualitatively new 
bounded solutions of 0. Therefore, violation of (Al) does not result in new types of wave 
solutions of Q satisfying (CI), (C2). 

Appearance or disappearance of u-fronts correspond to appearance or disappearance of the 
roots of the function h(u) which can occur with variation of the parameters c and a. The sim- 
plest case of the appearance of two or three roots corresponds to the fold or cusp bifurcations 
respectively |15| in the first equation of ([3]) . The simple conditions for the fold and cusp bifurca- 
tions show that, under variation of the boundary conditions (parameter a) and the wave speed 
(parameter c), appearance of traveling wave solutions of ([T]) is possible. 

Now we assume that (^43) does not hold, i.e., the functions gi(u) and h(u) have a coinciding 
root u. In this case system ([3]) has a line of non-isolated singularities in the phase plane (u,v). 
Each point of the form (u, v) is a non-isolated singular point; all the points on the line u = u 
are either simultaneously attracting or repelling in a transversal direction to this line [H] . If we 
assume that there exists a node (ui,v) of ([3]) such that u\ is a root of h(u), gi(ui) ^ 0, and there 
are no other roots of h(u) between u and u\ then, due to continuity arguments, there exists a 
family of bounded orbits of (Fig. [3|). To describe the traveling wave solutions corresponding 
to this family we define 

Definition 2. We shall say that model ([T]) possesses a family of free-boundary wave fronts in 
v-component if a) every v(z) — > v when z — > oo (z — > — oo); b) there exists an interval (^1,^2) 
such that for any v* £ (t>i, ^2) it is possible to find v-profile with the property v(z) — > v* when 
z — > —00 (z — > 00). 

Summarizing we obtain 

Proposition 3. The system (pQ) satisfying (CI), (C2), (^41), (A2) has a traveling wave solution 
such that u-profile is a front and v-profile is a free boundary front if and only if condition (A3) 
is violated and there is a node of system ([3]) such that there are no other singular points of ([3|) 
between this node and the line of non-isolated singular points. 
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The primary importance of such traveling wave solutions comes from the fact that for an 
arbitrary boundary condition (from a particular interval) for system (pQ) we can find a wave 
solution whose w-profile is a front. Note that violation of (A3) and simultaneous appearance of a 
free-boundary v component naturally occurs when the roots of h(u) are shifted under variation of 
c and a. It follows from ([3]) that bifurcation of w-profile occurs at (c*, a*) such that u = —a*/c* 
is a simple root of gi(u). 



2.3 An illustrative example 

Here we present a simple example to illustrate the theoretical analysis from the previous sections. 
We consider model (pQ) with 

f(u,v)= ( "~/ )(1 ~ u) , g(u,v) = -k(u-r)v(v-l), (7) 
v{v — 1) 

where I, k, r are non-negative parameters. The particular form of the functions f (u, v) , g(u, v) 
obviously satisfies (CI) and (C2). 
The wave system reads 

u = cu — k(u — l)(l — u)(u — r) + a, 
v' = —k(u — r)v(v — l)/c. 

The system ([8]) can have up to six singular points. For instance, if we fix the parameter values 
I = 0.2, r = 0.4, k = 1, a = —0.1, c = 0.3 then system ([8]) possesses six singular points; therefore, 
there are two orbit cells ensuring existence of traveling wave solutions of system (UJ . 

A phase portrait of (JH]) is shown in Fig. [H From Fig. H]it can be seen that, with the given 
parameter values, there exist two qualitatively different traveling wave solutions of the initial 
cross-diffusion system which correspond to two cases of Theorem [TJ 

Numerical solutions of system ([1]) with functions ([7|) and the given parameter values are 
shown in Fig. (the details of the numerical computations are presented in the Appendix). 

If we change the value of a to —0.12 then we obtain a family of free-boundary traveling wave 
solutions. 



3 Wave solutions of the separable models and some generaliza- 
tions 

3.1 General theory 

In this section we study models ([T|) which satisfy (CI). The rational function f2(v) = Z(v)/R(v) 
can be presented in the form 

f ( ) - Z ^ - Z(v)(v -Vl).-.(V- Vm) 

l2[V) ~ R(v) ~ R(v)(v-v 1 )...(v-v k ) , 
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Figure 4: The phase portrait of system (|8|). The parameters are I = 0.2, r = 0.4, k = 1, a 
-0.1, c = 0.3. 




Figure 5: Numerical solutions of system ([T]) with functions ©. The initial conditions are chosen 
to start the calculations from the orbit cells labelled / (panel (a), front-impulse solution) or 
II (panel (6), front-front solution) in Fig. UJ The solutions are shown for the time moments 
to = (bold curves) < t\ < ti < t% = 18 in equal time intervals 



where Z(v),R(v) do not have real roots; m ^ 0, k > 0; Vi ^ Vj for any The wave system 
has the form 

Ziv) 

cR(v) ' (9) 

v' = gi(u)g 2 {v)/c. 
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u' = cu + a + f l {u)g l {u)g 2 {v) 



By transformation of the independent variable 



dz , . 

* - -m- <10> 

which is smooth for any v except for v = Vi, i = l..k, system © becomes 

^ = c 2 R(v)(u + a/c) + fi(u)gi(u)g 2 (v)Z(v), 

I (») 

T = 9i(u)92(v)R{v). 
ay 

The roots of functions fi(u), gi(u), g2(v), Z(v), R(v) do not depend on parameters c and a, 
hence, we will suppose that the following conditions of non-degeneracy are fulfilled: 

(Bl) R(v) and g2(v) have no common roots; 

(B2) fi(u) and gi(u) have no common roots; 

(B3) fi(u),gi(u),g2(v),R(v) have no multiple roots. 

Coordinates of singular points (u, v) = (u, v) of (jlip can be found from one of the systems: 

-a/c = u, g 2 {v) = 0, (12) 
R(v) = 0, /i(u)=0, (13) 
R(v) = 0, gi(u) = 0, (14) 

or from combination of f)12[) - (|14[) . 

To infer possible types of the singular points of (llip we consider D = tr 2 J — 4 det J, where 
J is the Jacobian of (jlip evaluated at a singular point (u,v). 

If (u, v) is a solution of (TT2]) then 

tr J = J R(i))(c 2 + 5l ({i)^(i))), 
det J = c 2 5 i(u) J R(v) 2 ^(7)), 

D = J R(i)) 2 ( c 2 - 5l (n)^(l))) 2 ; 

if (n, v) is a solution of (fT3]) then 

tr J = 5l (n) 52 (t))(Z(t))/(^) + 

det J = ^(uJsaCfi)) 2 ^^)/!^)^^), 

D = {g l (u)g 2 {v)) 2 {Z{v)f' l {u)-R'{v)f; 

if (tt, u) is a solution of (|14p then 

detJ = 0, trJ = /i(% 2 («)^(%i(«)- 

Consequently we obtain that (u, v) is a saddle or a node for the cases corresponding to (fT2"|) 
and (TT3|) . and is a saddle, node, or saddle- node (see [H]) in the case (fT4"|) . Just as for the reduced 
separable models there are no singular points of (fTT|) of center or focus type. 
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Figure 6: Complex shape of wave profiles of model ©. The functions are f\(u) = u— 1, gx(u) = 
u, g2(v) = 1, R{v) = v(v — 1), Z(v) = 2v — 1. The left panel shows an impulse-front solution, 
the right panel shows a front-impulse solution 

Here we do not pursue the problem of classification of possible structures of the phase plane of 
([lip and are only concerned with new types of traveling wave solutions. Analyzing formulas for 
det J and tr J and applying arguments in the line with the proof of Proposition [2] and Theorem 
[Owe obtain 

Proposition 4. Let coordinates of singular points of (jllD satisfy (I12D and function gi (v) have 
two real neighboring roots v% and V2- If the function R(v) has an odd number of roots between 
v% and $2 and point (u = —a/c,v\) is a node, then (u = — a/c, v^) is a node as well (and vice 
versa) . 

Let coordinates of singular points of ([lip satisfy (JT3|) andu = u be a root of f\(u), R(v) have 
two real neighboring roots v% and 1)2 . If the function Z(v) has an odd number of roots between 
v\ and V2 and point (u,v\) is a node, then (n, V2) is a node as well (and vice versa). 

Due to the structure of system (jlip the phase plane is divided into horizontal strips, whose 
boundaries are given by v = v, where v is a root of R(v) or g2(v); all singular points of (jlip are 
situated on these boundaries. Bringing in the continuity arguments we obtain that under the 
conditions of Proposition H] there is a family of bounded orbits of (jlip which correspond to the 
traveling wave solution of ([T]) of an impulse-front type. 

It is worth noting that the structure of the phase plane inside a strip can be quite arbitrary, 
and we can only indicate asymptotic behavior of orbits in neighborhoods of singular points. As 
a result, under given boundary conditions (or, equivalently, fixed a) families of traveling wave 
solutions may have complex shapes (opposite to the examples presented in Fig. [T]). For instance, 
there can be non-monotonous fronts with humps or impulses which also have multiple humps 
and hollows. In general, we can only state that the form of impulses and fronts can be quite 
arbitrary, which is illustrated in Fig. [6j 

Under variation of parameters c and a it is possible that function g\ (u) has a root u = —a*/c* 
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for some values of the parameters; in this case system (lllh has a line of non- isolated singular 
points u = u and the analysis in this situation is similar to the analysis which led to Proposition 
El a family of free-boundary fronts appears in v component. 

Now let us relax the condition (B2); we assume that there exists such u that f\{u) = and 
g\(u) = 0. We can always find values of c and a such that u = —a/c. In this case the phase plane 
of (jlip has a line u = u of non-isolated singular points. After the change of the independent 
variable dr = dy/{u — u) the resulting system still possesses singular point of the form (u,v), 
where v is a root of R(v). If this point is a node then, applying continuity arguments, we obtain 
that there exist a family of orbits of system (|lip such that some solutions from a neighborhood 
of (u,v) tend to this point if y — > oo (or y — > — oo) and tend to point of the form (u, v*) if 
y — * — oo (or y — > oo), where v* is an arbitrary constant from some interval. These solutions 
correspond to traveling wave solution of ([T]) such that -u-profile is an impulse and v component 
is a free-boundary front (see Fig. [8|). 

Summarizing we obtain 

Theorem 2. The system ((TJ) satisfying (CI) and (B1)-(B3) can only possess traveling wave 
solutions of the following kinds: 

i. front-front solutions; 

ii. front-impulse solutions; 

iii. impulse-front solutions; 

iv. Under variation of parameters c and a it is possible to have wave solutions where u 
component is a front, v component is a free-boundary front; 

v. Under the additional condition that (B2) does not hold it is possible to have wave solutions 
with u component is an impulse and v component is a free-boundary front. 



3.2 The phase plane analysis of the Keller-Segel model 

The classical Keller-Segel model has the form ([T]) with f{u,v) = 5u/v, g(u,v) = —ku, where 
5, k > (see [3j). In our terminology this model falls in the class of the separable models; the 
wave system reads 

u 2 

u = cu — 5k h a, , 1 

cv (15) 

v' = —ku/c, 

which, with the help of transformation ()10p . can be reduced to the form (|11|) : 

dU n 2 

— = c uv — oku + acv, 

t „ < 16 > 

— = —kuv. 
dy 

If a ^ then system (|16|) has the only singular point (it, v) = (0, 0). In this case the separable 
model cannot possess traveling wave solutions (Theorem [2]). Hence we put a = 0. Note that 
the last requirement is necessary if one supposes that U(x,t) should be finite. 
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For a = system (I16D has a line of non-isolated singular points (0, v), and there is also 
additional degeneracy at the point (0,0) (case v. in Theorem ([2]))- This can be seen applying 
the second transformation of the independent variable dr = udy, which leads to the system 




(17) 



for which the origin is a topological node with the eigenvalues Ai = — k5 and A2 = —k. Thus 
there exists a family of traveling wave solutions whose u-profile is an impulse and u-profile is a 
free-boundary front. 

In Fig. [7] we show how the parametrization of the phase curves of the wave system change 
after the transformations of the independent variables. This picture can also serve as an illus- 
tration to assertion v. of Theorem [2j 

Due to biological interpretation of the Keller-Segel model it is necessary to have U(x, t) ^ 
and V(x,t) ^ 0. Using the fact that an eigenvector (1,(<5 — l)k/c 2 ) corresponds to Ai and 
(1,0) corresponds to A2 it is straightforward to see that to ensure the existence of non- negative 
travelling way solutions we should have 5 > 1. 

Numerical solutions of the Keller-Segel model are given in Fig. [8b. Originally, the Keller- 
Segel model was suggested to describe movement of bands of E. coli which were observed to 
travel at a constant speed when the bacteria are placed in one end of a capillary tube containing 
oxygen and an energy source [3j. In Fig. [8b it can be seen that bacteria (u(x,t)) seek an 
optimal environment: the bacteria avoid low concentrations and move preferentially toward 
higher concentrations of some critical substrate (v(x, t)). The stability of the traveling solutions 
found was studied analytically in [I9j [21] . 

3.3 Impulse-impulse solutions 

In the preceding sections we studied the systems (P) where the functions f(u,v),g(u,v) can be 
represented as a product of functions that only depend on one variable. The next natural step 




Figure 7: The phase planes of systems (fTTl) (a), (fl"6l) (6) and (fT5j) (c) 
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U -60 -40 -20 x 

Figure 8: (a) The phase plane of system (fT5|k the parameters are a = 0, k = 1, 5 = 4, c = 0.5. (6) 
Numerical solutions of the Keller-Segel system with the parameters given in (a). The solutions 
are shown for the time moments to = (bold curves) < t\ < t% < t% = 30 in equal time intervals 



is to assume that these functions depend on affine expressions au + bv + c, where a, b, c are 
not equal to zero simultaneously. Here we present an explicit example of such a system. The 
example is motivated by appearance of a particular type of traveling wave solutions, which is 
absent in the separable models. 
We suppose that 

Su 

f{u,v) = ——, 5>0,P>0, 

(3u + v (18) 

g(u, v) = —ku + rv, k, r ^ 0. 

If P, r = then we have the Keller-Segel model studied in Section 3.2. 
The wave system for (|TJ) with the functions given by (|18p reads 

, 5u(—ku + rv) 

U = CU-\ 7- r , 

c(pu + v) (19) 
v' = (—ku + rv)/c, 

where we put parameter a equal zero. 

After the change of the independent variable dz / (c(Pu + v)) = dy system (|19|) takes the form 

— = c 2 u(Pu + v) + 5u(—ku + rv), 

I < 2 »> 

— = (—ku + rv)(Pu + v). 
dy 

If r = then system (|20p has a line on non-isolated singular points u = 0; if r ^ then 
(u,v) = (0,0) is an isolated non- hyperbolic singular point of (f20l) (i.e., both eigenvalues of the 
Jacobian evaluated at this point are zero). 
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Figure 9: (a) The phase plane of system (fT9j) : the parameters are k = 1, j3 = 1, 8 = 2, c = 
0.43, r = 0.1. (b) Numerical solutions of system ([1]) with the functions given by (|18p and the 
parameter values as in (o). The solutions are shown for the time moments to = (bold curves) < 
t\ < t<i < t% < £4 < £5 = 20 in equal time intervals 



First we consider the case r = 0. After yet another transformation dr = dy/u we obtain 

(c 2 /3 - 5k)u + c 2 v, 

(21) 

-k((5u + v). 

Thus for 5 > 1 and 0^/3^ k{^/~8— l) 2 /c 2 the origin is a node for system (|2ip which implies that 
in the initial system (|20p there exist a family of bounded orbits which tent to (0, 0) for y — » 00 
or y — > — 00. This family corresponds to a family of traveling wave solution of the system ([T|) 
where n-profile is an impulse and u-profile is a free-boundary front. The picture is topologically 
equivalent to the phase portrait shown in Fig. [Hh. 

For r / the wave system (]20p has singular point (0, 0) possessing two elliptic sectors in its 
neighborhood (see Fig. [^)- The proof of existence of the elliptic sectors can easily be conducted 
with the methods given in [24j . Asymptotics of homoclinics composing the elliptic sector are 
u = (trivial) and v = K + u, where K + is the biggest root of the equation 

K 2 (c 2 + r(S - 1)) + K(pc 2 - k(5 - 1) - f3r) + 0k = 0. 

The family of homoclinics in the phase plane (it, v) correspond to the family of wave impulses 
for the system ([1]) (see Fig. [lb and Fig. [9b>). To our knowledge such kind of solutions (infinitely 
many traveling wave solutions of impulse-impulse type with the fixed values of the parameters) 
was not previously described in the literature. 

The results of the numerical computations (Fig. Efc>) indicate that the family of traveling 
impulses is clearly non-stationary, since its amplitude decays visibly in time. Which is important, 
however, is that it is possible to observe moving impulses at least at a finite time interval. 



du 

J? 
dv 
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4 Conclusions 



In this paper we described all possible traveling wave solutions {u(z), v(z), z = x + ct} of 
the cross-diffusion two-component PDE model (pQ) satisfying (CI), where the cross-diffusion 
coefficient may depend on both variables and possess singularities. Such kind of models is 
widely used in modeling populations that can chemotactically react to an immovable signal 
(attractant) (e.g., [BE] and references therein). The study of traveling wave solutions of model 
(PQ) was carried out by qualitative and bifurcation analysis of the phase portraits of the wave 
system ([2]) that depends on parameters c and a. Here c is a speed of wave propagation and 
a characterizes boundary conditions under given c. Any traveling wave solution of (pQ) with 
given boundary conditions corresponds to a solution of wave system ([2]) with specific values 
of parameters c and a; the converse is also true. Therefore, instead of trying to construct a 
traveling wave solution with the given boundary conditions, we study the set of all possible 
bounded solutions of the wave system considering c and a as its parameters. This approach 
allows identifying all boundary conditions for which model (pQ) possesses traveling wave solutions. 

The main attention is paid to the so-called separable model, i.e., to model ([1]) that satisfies 
(CI); it is worth noting that in this case the functions f(u,v) and g(u,v) in (pQ) are products of 
factors that depend on a single variable. We showed that for some fixed values of parameters 
c and a the solutions of the wave system compose two-parameter family (for explicit formulas 
in a simple particular case see ([5])). One result is that {u(z), w(z)}-profiles of the wave solution 
of the reduced separable models (pQ) that satisfies (CI), (C2) can be only front-front or front- 
impulse; for more general case of the separable models we can additionally have impulse-front 
solutions. For some special relations between c, a and the model parameters model (pQ) can 
have wave solutions whose v -profiles are free-boundary fronts, i.e., v(z) tends to an arbitrary 
constant from some interval at z — > oo or at z — ► — oo. Note that traveling wave solutions of the 
well-known Keller-Segel model [3] as well as of its generalization (Section 3.3) have impulse-front 
profiles with a free-boundary front (Sections 3.2, 3.3). 

We also considered a natural extension of the separable models; namely, we gave an example 
of model (pQ) where /, g are products of factors that depend on affine expression of both variables 
(see (|18p ). This model can be considered as a generalization of the Keller-Segel model because it 
has two additional parameters and turns into the Keller-Segel model if both of these parameters 
are zero. If only one of the parameters vanishes, the model has a family of "Keller-Segel" -type 
solutions, i.e., u-profile is an impulse and u-profile is a frond with a free boundary. Importantly, 
in some parameter domains this model possesses a two-parameter family of impulse-impulse 
solutions (Fig. [9]). To the best of our knowledge, such type of traveling wave solutions was not 
previously described in the literature: depending on the initial conditions traveling impulses can 
have quite a different form for the fixed values of the system parameters. Taking into account 
the fact that such solutions are absent in the separable models, we can consider model ([1]) with 
functions (|18|) as the simplest model possessing this type of traveling wave solutions. 

Rearrangements of traveling wave solutions of PDE model (pQ), which occur with changes of 
the wave propagation velocity and the boundary conditions, correspond to bifurcations of its 
ODE wave system. In particular, appearance/disappearance of front-profiles with variation of 
parameters c and a correspond to the fold or cusp bifurcations in the wave system; rearrange- 
ment of a front to an impulse can be accompanied by appearance/disappearance of non-isolated 
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singular points in the phase plane of the corresponding wave system (see Section 2.3). Existence 
of non-isolated singular points in the wave system may result in the existence of free-boundary 
fronts in model (fT]). For instance, this is the case for the Keller-Segel model. 

We emphasize that the separable model ([1]) , when the values of parameters c and a fixed in 
the wave system, possesses, in general, two-parameter family of traveling wave solutions. There 
are infinitely many bounded orbits of ([2]) that correspond to traveling wave solutions of (pQ) (see 
Figures [5J O El ©• It is of particular interest that in all numerical solutions of ([T]) that we 
conducted it is possible to observe traveling waves. We did not discuss the issue of stability of 
the traveling wave solutions found, but we note that it is usually true that unstable solutions 
cannot be produced in numerical calculations. It is tempting to put forward a hypothesis that 
the presence of additional degrees of freedom (two free parameters) is the reason of producing 
traveling waves in numerical computations. This important question can be a subject of future 
research. 



5 Appendix 

We did numerical simulations of system ([I]) for x £ [-L, L], where L varied in different numerical 
experiments. We used no-flux boundary conditions for the spatial variable. Inasmuch as we 
wanted to study the behavior of the traveling wave solutions in an infinite space we chose such 
space interval so that to avoid the influence of boundaries. 

We used an explicit difference scheme. The approximation of the taxis term is an "up- 
wind" explicit scheme [25] which is frequently used for cross-diffusion systems (e.g., [26]). More 
precisely, 



v^ 1 =vl + (At)g(uivj), i = 2,...,N-l, 



where for the positive taxis (pursuit) (i.e., f(u,v) < 0), 

A = f(ul,vj) if Vi+i > Vi, 

A = / («*+!, vj +1 ) if v i+1 < Vi, 

B = f( ut i-li v l-l) if v i > v i-!, 



B = f(u\,v t i) iivi<v. 



i-l- 



For the negative taxis (invasion): 

A = f{u\ 1 v\) iiv i+l <v,i 
A = f{u\ +1 ,v\ +1 ) if v i+ i ^ n 

B = /(i4-i.«f-i) if«i<«j_i 
B = f(u\,v\) Hvi^Vi-i. 

We used At = 0.001, Ax = 0.1. For the boundary conditions: 

,.t+i _ ,.t „*+i _ „t 

7,* + 1 - if V t+1 - 7,* 

^i —v 2 , v N — Ujv_i- 
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For the initial conditions we used numerical solutions of the corresponding wave systems. 
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